Absolute timing of the Crab pulsar with the INTEGRAL/SPI telescope 



S. Molkov, E. Jourdain and J. P. Roques 
molkov@iki . rssi . ru 

CESR - Universite de Toulouse (UPS), CNRS (UMR 5187), 9 Av. du Colonel Roche, 31028 Toulouse Cedex 4, France 



ABSTRACT 

We have investigated the pulse shape evolution of the Crab pulsar emission in the hard X-ray domain 
of the electromagnetic spectrum. In particular, we have studied the alignment of the Crab pulsar phase 
profiles measured in the hard X-rays and in other wavebands. To obtain the hard X-ray pulse profiles, we 
have used six year (2003 - 2009, with a total exposure of about 4 Ms) of publicly available data of the SPI 
telescope on-board of the INTEGRAL observatory, folded with the pulsar time solution derived from the 
Jodrell Bank Crab Pulsar Monthly Ephemeris (http://www.jb.man.ac.uk). We found that the main pulse 
in the hard X-ray 20 - 1 00 ke V energy band is leading the radio one by 8 . 1 8 + 0.46 milliperiods in phase, 
or 275 + 15 /us in time. Quoted errors represent only statistical uncertainties. Our systematic error is 
estimated to be ~ 4-Q/is and is mainly caused by the radio measurement uncertainties. In hard X-rays, the 
average distance between the main pulse and interpulse on the phase plane is 0.3989 ±0.0009. To compare 
our findings in hard X-rays with the soft 2-20 keV X-ray band, we have used data of quasi-simultaneous 
Crab observations with the PCA monitor on-board the Rossi X-Ray Timing Explorer (RXTE) mission. 
The time lag and the pulses separation values measured in the 3-20 keV band are 0.00933 + 0.00016 
(corresponding to 310 + 6 /us) and 0.40016 + 0.00028 parts of the cycle, respectively. While the pulse 
separation values measured in soft X-ray s and hard X-ray s agree, the time lags are statistically different. 
Additional analysis show that the delay between the radio and X-ray signals varies with energy in the 2 
- 300 keV energy range. We explain such a behaviour as due to the superposition of two independent 
components responsible for the Crab pulsed emission in this energy band. 

Subject headings: stars:neutron-pulsars: general-pulsars: individual: PSRB0531+21 - X-rays:star 



1. Introduction 

The Crab pulsar (PSR B0531+21) is the best stud- 
ied isolated pulsar. The pulsed e mission was dis- 
covered long ag o and in the X-ray s dFritz et al. 1969: 
Bradt et al. 1969t) and in the y-rays dKurfess 197 lb do- 
mains, and its pulse morphology has been studied in 
the full range of the electromagnetic spectrum. In all 
energies, the pulse profile has two prominent features, 
the main pulse (or the first peak, PI) and the interpulse 
(or the second peak, P2). The relative intensities of 
these peaks depend on the energy band. The second 
peak dominates in the ~ 200 - 1200 keV energy band. 



Based on observations with INTEGRAL, an ESA project with 
instruments and science data centre funded by ESA member states 
(especially the PI countries: Denmark, France, Germany, Italy, 
Switzerland, Spain), Czech Republic and Poland, and with the par- 
ticipation of Russia and the USA. 



By all appearances, the distance between the peaks 
on the phase plane is almost constant in time, slightly 
varying aro und the value Ai/r - 0.40(0) depending on 
energy (see iThompson et al. 19771: I Wills etal. 19821: 
White et al 1985HNolan et al. 1993: Masnou et al. 1994 



Moffe tt and Hankins 1996UPravdo et al. 1 9971: iKuiper et al. 200lt 



Brandt et al. 2003l:lRots et al. 2004tlthe MAGIC collaboration 2 008) 
For a long time it has been assumed that both peaks 
are perfectly lined up in phase over the whole en- 
ergy range. This assumption has been disputed for 
the first time in the work of Masnou et al. (1994). 
Based on the data of the FIGARO II telescope (bal- 
loon experiment), the authors found that the first peak 
in the 0.15-4 MeV energy band is leading the radio 
main pulse by ~ 400 + 150 fis. Later, the misalign- 
ment in phase of the main radio pulse and the main 
pulse in shorter wavelengths has been confirmed by 
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several instruments. No absolute agreement exists in 
the value of the radio delay measured by different in- 
struments even for the same energy band, though they 
are close to each other especially if one takes into ac- 
count not only statistical errors but also the possible 
systematic uncertainties. The most recent measure- 
ments of the X ,y-rays to radio l ag are: RXTE/PCA — 
344 ± 40 fj.s dRots et al. 20041). JEM- X/INTEGRAL 
— 300 + 67 us dBrandt et al. 20031) . ISGRI/ and 
SPI/INTEGRAL — 285 ± 12 us and 265 + 23 us, re - 
spectively (statistical errors only, [ Kuiper et al. 2003), 
EGRET/CGRO — 241 + 29 p.s (IKuiper et al. 20031) . 
The INTEGRAL results mentioned above are based 
on observations covering only several days, thus these 
data can be folded with the single ephemeris record, 
while different ephemerides were used for RXTE 
and CGRO observations. That means that RXTE 
and CGRO results should be statistically more sig- 
nificant since they are less af fected by uncert ainties 
in the radio ephemerides (see iRots et al. 20041) . Re- 
cently, the opti cal-radio delay has also been confirmed, 
255 + 21 fis dOosterbroek et al. 20081) . 

In this paper we present the timing analysis of the 
Crab pulsar in the 2 - 300 keV energy range with the 
SPI/INTEGRAL telescope and the RXTE instruments. 

2. Observations and data reduction 
2.1. INTEGRAL 

The International Gamma-Ray Astrop hysics Labo- 
ratory (INTEGRAL; IWinkler et al. 2003 was launched 
on Proton/Block-DM on October 17, 2002 into a 
geosync hronous highly ecc entric orbit with high 
perigee dEismont et al. 20031) . The scientific pay- 
load of INTEGRA L includes four te lescopes: the 
spectr ometer SPI ( Vedrenne et al. 20031) . the imager 
IBI S dUbertini et al 20031) . the X-rav monitor JEM- 
X dLund et al. 20031) a nd the optical monitor OMC 
dMas-Hesse et al. 20031) . 



2.7.7. The SPI telescope characteristics 

In this paper, we focus on analysis of the SPI data. 
For detailed description of the instrument calibrations 
and performance, see Vedrenne et al. (2003), Attie et 
al. (2003) and Roques et al. (2003). Below we give 
only a brief review of the key characteristics of the in- 
strument and software relevant for this work. 

SPI consists of 19 high purity germanium detectors 
(GeD) packed into a hexagonal array (see Figure Q]). 



Fig. 1 . — The shadowgram of the detectors plane of the 
SPI telescope in the simple case when only one source 
is located in the SPI FOV. The"black" colour corre- 
sponds to the area illuminated by the source ; "Dark 
grey" to the area shadowed by the mask; "Light grey - 
to the area closed for the source by the collimator. 



The combination of a good sensitivity to the contin- 
uum and line emission in the energy band 20 - 8000 
keV, provided by a large geometrical detector's plane 
area (~ 500 cm 2 ) and cryogenic system, and a good 
timing resolution (102.4 fisec) gives us an opportunity 
to study even very fast X/y-rays pulsars. 

There are two main types of SPI events: single 
events (registered by only one detector) and multiple 
events (scattered photon detected by two and more de- 
tectors). Though the pulsating signal is clearly de- 
tected in both types of events, we did not use the mul- 
tiple events because of difficulties in the extraction of 
the spatial information, worse timing resolution and 
the lack of a tested energy response matrix. It should 
be noted that this kind of data was successfully used 
by Dean et al. (2008) for measuring Crab polarisation. 

2.1.2. The SPI clock and time connection 

SPI has its own internal clock, a 20MHz oscillator 
that generates time tag signal every 2048 periods, i.e. 
102.4 jjs - the time resolution of the instrument. The 
SPI clock is synchronised with the on-board clock at 
every 125 ms (the 8Hz on-board cycle, OBT%hz) by a 
reseting of the counter associated to the SPI oscillator. 



2 



The SPI time tags are counted from the beginning of 
the 8Hz cycle and divide the 125 ms interval into 1019 
intervals of 102.4 micro seconds duration plus one in- 
terval with a duration of 70 micro seconds. In fact, 
the base frequencies of both on-board and SPI clocks 
are slightly varying with time (e.g. due to variation 
of the crystals temperature) that leads to variation in 
duration of SPI time tag intervals. We have estimated 
this effect to be very small (< 5 [ts for any statistically 
significant analysis), neglected it and used the nominal 
values of the frequencies. The on-board time (OBT) 
of the current 8Hz cycle is the value of the datation of 
the beginning of the cycle (SPI User Manual - Issue 
5.2 - SEP 2002). By convention, the On Board Time 
(OBT) of any event is the time of the leading edge of 
the interval where the event has been detected, that is : 



OBT = OBT mz +N„*D/2 



-20 



(1) 



where N„ is the number of the SPI time tag interval 
(0..1019) and D = 102.4 x 10 ~ 6 s is the length of the 
interval. The divisor 2~ 20 is introduced to convert time 
expressed in seconds to on-board time units and re- 
flects the conventional time resolution of the INTE- 
GRAL clock (the real clock accuracy is 2~ 19 seconds). 
For the conversion of the On Board Time to Coor- 
dinated Universal Time (UTC) and for the barycen- 
tric correction we have used routines from the stan- 
dard Off-line Science Analysis software package ver- 
sion 7.0 developed at INTE GRAL Science Data Cen- 
tre (ICourvoisier et al. 20031) and the time correlation 
files provided with the auxiliary data for the INTE- 
GRAL data archive generation number 2. The equa- 
tion (1) and the time transformation routines ensure 
accuracy of returning Universal Time of the order of 
100 ids. For precise timing analysis, a time correction 
must be added to OBT before any conversion: 

AT = AT SPI + AT sal + AT rev2 , (2) 

where AT$pi = D/2 = 51.2 fis — the mean system- 
atic shift due to the fact that the arrival time of events 
is defined as the time of the leading edge of the time 
tag interval, while the actual arrival times are normally 
distributed inside the time tag interval; 
AT sat = 83 fj.s - the delay betwee n the on-board time 
and SPI time (ground calibration, lAlenia Spazio 2 002): 
ATrevi = -47 fis - this shift is common for all time 
correlation files from the archive generation 2 (will 
disappear in the new generation of archive). 
Note, that the SPI instrumental delay given in the Ta- 
ble 4 of Walter et al. (2003) is the sum of the first two 
terms (134 fis) of the expression (2). 



In our analysis, to avoid an additional discretiza- 
tion in time series we decided to define the time of any 
registered event not as the time of the leading edge of 
the appropriate time tag interval, but as a linearly ran- 
domised time inside this time tag interval, and the final 
expression for OBT is the following: 

OBT = OBT %hz + ((N„ + rand[0, l])*D + AT) /2~ 20 , 

(3) 

in this case ATspi = and therefore AT = +36 x 
10~ 6 s. 



Table 1: List of observations used in this analysis. 



Rev. 


Observing Period 


Exp." 


Target * 


N° 


UTC 


ks 




YEAR 2003 


0043 


19.177 


- 21.783 Feb. 


164 


CRAB 


0044 


22.163 


- 24.772 Feb. 


177 


CRAB 


0045 


25.167 


- 27.699 Feb. 


162 


CRAB 


0102 


14.623 


- 17.218 Aug. 


91 


CRAB 


0103 


17.610 


- 17.845 Aug. 


22 


CRAB 


0124 


19.430 


- 22.031 Oct. 


195 


IC443 


0125 


22.422 


- 25.023 Oct. 


196 


IC443 


0126 


25.414 


- 28.015 Oct. 


110 


IC443 


YEAR 2004 


0170 


5.122 - 


7.302 Mar. 


119 


CRAB 


0182 


9.987 - 


12.504 Apr. 


197 


IC443 


0184 


15.976 


- 18.566 Apr. 


222 


IC443 


0239 


27.435 


- 30.036 Sep. 


187 


CRAB 


0247 


21.363 


- 23.964 Oct. 


223 


IC433 


YEAR 2005 


0300 


28.910 


- 31.504 Mar. 


188 


CRAB 


0352 


31.463 


- 2.890 Aug/Sep. 


200 


A0535 


0365 


11.096 


- 11.938 Oct. 


59 


CRAB 


YEAR 2006 


0422 


28.729 


- 31.198 Mar. 


190 


CRAB 


0464 


1.252 - 


3.011 Aug. 


60 


Taurus 


0483 


28.711 


- 29.627 Sep. 


69 


CRAB 


YEAR 2007 


0541 


19.537 


- 22.128 Mar. 


213 


CRAB 


0605 


27.016 


- 28.927 Sep. 


154 


CRAB 


YEAR 2008 


0665 


24.473 


- 27.033 Mar. 


194 


CRAB 


0666 


27.463 


- 30.022 Mar. 


204 


CRAB 


0727 


25.940 


- 28.488 Sep. 


197 


CRAB 


0728 


28.932 


- 1.499 Sep./Oct. 


184 


CRAB 


YEAR 2009 


0774 


13.538 


-16.145 Feb. 


196 


CRAB 



" - this value represents the total exposure of selected data 
b - for the complete description of the observations see the 
ISOC site http://www.sciops.esa.int 
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2.1.3. Data selection 

In our analysis we have used publicly available data 
of all observations where the source was in the FOV of 
the SPI telescope. The exposure of selected data totals 
up approximately 4 Ms. 

SPI is a telescope with a coded mask aperture and 
most INTEGRAL observations are organised as a set 
of snapshots of the sky around a target (pointings or 
Science Windows - continuous observations pointed 
on a given direction in the sky). It means that the in- 
strument effective area for the chosen target is chang- 
ing from pointing to pointing. This effective area is 
mainly determined by geometrical area of the non- 
shadowed part of the detector plane (see Figure[T|i. 

To extract the Crab pu lsed signal we hav e used the 
epoch folding technique (ILeahyetal. 19831) . Any set 
of observations can be represented as a set of the whole 
detector plane count rates and in these terms, the to- 
tal folded curve is the direct sum of the folded coun- 
trates of the individual pointings. To reach the best 
result we need to find an optimal series of the INTE- 
GRAL pointings for which the signal to noise ratio will 
be the highest one. In this regard, several parameters 
characterise each pointing: the exposure of the point- 
ing, Tf, the background conditions - the instrumental 
background countrate plus the sum of the countrates 
from other sources in the FOV, C b ; the illumination 
fraction of the detector plane, a,-, corresponding to the 
source direction (could include not only geometrical 
factor, varies from to ~ 0.6); and the mean countrate 
of the pulsating part of the source emission, C p (for 
completeness, below we are introducing also the term 
C dc - the mean unpulsating countrate of the pulsar, for 
the Crab pulsar C dc = 0). Using the terminology in- 
troduced above and assuming that all variances follow 
the Poisson statistic, the signal to noise ratio for the 
sequence of M pointings can be expressed as follows: 



\N)m 



M 



M i M 

E(cf + a,(cf + cf e ))r« 



(4) 



The optimal set of K pointings chosen from the ini- 
tial set of M pointings is that for which the value of 
{S IN)k calculated using the equation (4) reaches its 
maximum. In the case of a large number M the exhaus- 
tive search will take infinite time (even for M — 100 
the number of combinations exceeds 10 30 , in our case 




15000 30000 
K 



Fig. 2. — The left panels show the distribution of the 
number of SCWs (a) or the number of individual de- 
tectors (b) versus their illumination fraction for our 
data set. The function F(K) (see the equation (6)) is 
plotted on the rights panels for two cases: a) the whole 
detector plane; b) individual detectors. Shaded areas 
correspond to the data that should be excluded from 
the analysis. The function F in the case (b) is renor- 
malised to be in the same units as in the case (a). 



M > 1000). Several simplifications have been investi- 
gated. 

Considering that the INTEGRAL pointings (or Sci- 
ence Windows, SCWs) cover generally ~ 2 - 5 ks 
time intervals, we can treat the exposure T, of in- 
dividual member of the data series as constant, i.e. 
Ti =a const - T, i e [1,M]. Moreover, in the SPI 
telescope, the background dominates the useful signal, 
i.e. (C p + C. ) «: C b is true for any point source ex- 
cluding the brightest events like Gamma-Ray Bursts or 
short intense bursts from Soft Gamma-Ray Repeaters 
and Anomalous X-ray Pulsars. Further, the amplitude 
of the variation of the SPI background does not ex- 
ceed ~ 50% of its mean level (excluding periods of 
Solar flares). So, we can treat C b - const - Cb for 
i G [1,M]. We are assuming also that the mean value 
of the pulsar "Pulsed" component is constant in time 
Cf - const - C p . Now, taking into account assump- 
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tions listed above, we can modify the equation (4) as: 



\N/m \ 



M 
1=1 



(5) 



Cb Vm 

In equation (5) we do not reduce the number of com- 
binations in comparison with the equation (4) but now 
it is easy to see that the procedure of searching for the 
optimal set is equivalent to searching for the maximum 
value of the following discrete function: 



F(K) = 



K 

i=i 

"vr 



Ke [l,M] 



(6) 



where B = \J3i] is the back ordered A = [or,-] set. 

The speculations presented above could be easily 
extended to the case when we do not treat the whole 
detector plane (hereinafter referred to as the case I) but 
each detector separately (case II). In the latter case, a, 
is the illumination fraction of an individual detector 
and varies in the range 0-1, and M is the number of 
pointings times the number of the individual detectors 
(19 for SPI). 

We have obtained the solution for our set of INTE- 
GRAL observations for both cases as illustrated in Fig- 
ure [2] To reach the maximum of the signal to noise 
ratio in case I, we should exclude from the analysis the 
pointings with an illumination fraction below 21%. In 
case II, we should use only those detectors that have an 
illumination fraction above 38%. Figure|2l shows that 
using individual detectors (case II) we obtain a 30% 
improvement of the signal detection significance. In 
this paper dedicated to the timing analysis, we have 
implemented case II. 

2.2. Jodrell Bank Crab Pulsar Monthly Ephemeris 

For the folding of the Crab pulsar lightcurves, we 
have used the time solution derived from Jodrell Bank 
Crab Pulsar Monthly Ephemeris dLyne et al. 1993 ) 
and the corresponding Crab Pulsar coordinates. The 
database is available through the World Wide Web 
(http://www.jb.man.ac.uk/pulsar/crab.html) and con- 
tains the dispersion-corrected time of arrival of the 
centre of the main pulse (in TDB time system), the fre- 
quency and its first derivative and the range of validity. 
From this database we extracted the radio ephemerides 
covering the periods of the INTEGRAL observations. 
For each radio ephemeris record we took its two neigh- 
bours and calculated the second derivative of the fre- 
quencies so that the phases and frequencies given at 




200 600 1000 

Science Window number 



Fig. 3. — Science window by Science window folded 
curves in the 20-100 keV energy band. The folding 
procedure is based on ephemerides from the Table 2. 
No A*!* correction is applied (see text). 



the edge of the validity intervals and those deduced by 
extrapolation are consistent with each other better than 
0.0001 in phase and 10~ 7 in frequency. The resulting 
ephemerides for the INTEGRAL observing periods 
are given in Table 2. 

The main pulse arrival time in the monthly ephemeris 
is determined with an error around 60 fis, that includes 
the uncertainty in the delay due to interstellar scatter- 
ing ~ 20 fis (owing to the dispersion measure uncer- 
tainty ADM ~ 0.005 pc/cm 3 ) as well as those arising 
from unknown instrumental effects ~ 40 fis (see e.g. 
Rots et al. 2004 and references there). While the first 
part can be treated as a statistical error that follows 
the Poisson statistic and decreases with the number of 
independent measurements, the second part should be 
treated as a systematic error that is always present in 
the measured values. 

2.3. RXTE 

The PCA instrument onboard the RXTE (Rossi X- 
ray Timing Explorer) orbiting X-ray observatory con- 
sists of five identical proportional counters with a total 
area o f 6500 cm 2 , oper ating in the 2-60 keV energy 
range (IBradt et al. 1993b . The accuracy of the RXTE 
clock in absolute time for our observing time interval 
(2003-2009 yy.) is better than 2 ps (see Rots et al. 
2004 and references therein). Because of its large area 
and excellent time resolution and time accuracy, the 
instrument is sensitive enough to reconstruct a signif- 
icant 400 bins phase curve of the Crab pulsar using 
an exposure of the order of 1 ks. We used Crab PCA 
observations that coincide in time (within two weeks) 
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Fig. 4. — Crab phase histogram in the 20-100 keV en- 
ergy band in absolute phase with the phase resolution 
of 0.0025. "1" corresponds to the phase of the main 
radio pulse. Inner panels magnify the main pulse and 
the interpulse peaks. Dashed vertical lines bound the 
intervals used for the fit procedure. The exposure is of 
the order of 100 ks. 

with any of our INTEGRAL observation and contain 
data in the Generic event mode format with time reso- 
lution better than 250 fis. For the fine clock correction 
and barycenter correction we used faxbary script from 
the FTOOLS package that calls the axBary code (see 
e.g. the RXTE Guest Observer Facility). For the fold- 
ing procedure we use the same routine and the same 
ephemerides as for SPI/INTEGRAL. 

3. SPI/INTEGRAL analysis and results 

To perform epoch folding analysis we ascribe to 
each detected photon the phase using the appro- 
priate ephemeris from the Table 2 and the following 
formula: 



V t = V + f{t-t ) + \{t-t f dl ■ • > 3 



where to, /, -jt is the radio ephemerides valid 
for the moment "t" and = (we want to work in 
absolute phase, i.e. the main radio pulse is at phase 
0.0). Then we can plot the phase values producing 



light curves with the requested resolution (here we use 
400 bins per cycle). The middle of the zero bin is cor- 
responding to the phase 0.0. 

As a first step, to check our SPI data set for the 
presence of unknown "glitches" from the Crab pulsar 
or some instrumental artifacts (e.g. inaccuracy in the 
on-board to Universal time conversion procedure) we 
folded separately each SCW lightcurve in the broad 
20 - 100 keV energy band. The result of the dynami- 
cal folding is presented on Figure[3] where we see that 
the shape of the Crab phase histograms and absolute 
phases are very stable (the two peaks of the pulse are 
good tracers). 

In order to determine the phase of the hard X-ray 
main pulse and interpulse more precisely and to study 
possible variations of these values in time we summed 
up the SCW folded curves for each revolution. The SPI 
pulse profile for the Crab pulsar in the 20 - 100 keV 
energy band for the exposure ~ 100 ks is presented on 
Figure [4] To define the phase of the pulses, we fitted 
the data for each revolution with a composite model: a 
Gaussian function plus a constant background, in the 
phase intervals 0.98 - 1 .01 and 1 .38 - 1 .41 for the main 
pulse and interpulse, respectively, and adopted the fit- 
ted position of the Gaussian centroid as the appropriate 
phase of the corresponding pulse. To make sure that 
the fit results are model independent, we fitted the data 
for the main pulse with two other models: Lorentzian 
— used in Rots et al. (2004) (due to the lack of statis- 
tic for the one revolution timescale we could not apply 
the complete procedure of the peak-finding described 
in this paper) and Lorentzian plus constant — used in 
Kuiper et al. (2003). We found that all three mod- 
els yield similar values, with dispersion not exceeding 
0.0003 period (< 10 fis). 

The absolute phases of the main pulse and inter- 
pulse (with respect to the radio main pulse) together 
with the phase difference between them, as a func- 
tion of the observation number, are shown in Figure 
The distributions of *Ppi and *Pp2 values are well 
consistent with normal distributions Q¥p2 - T^i is not 
a measured value but a combination of two indepen- 
dent functions). The mean value of the hard X-ray 
main pulse phase is 0.99182 + 0.00046, thus, the pulse 
leads the radio main pulse by 8.18 ± 0.46 milliperi- 
ods or 275 ± 15 fis. Such quoted error does not in- 
clude 40 lis coming from the uncertainty of the radio 
ephemeris. The average value of the two pulses sepa- 
ration is 0.3989 ± 0.0009 parts of the cycle. Note, that 
this relative value is independent of any uncertainty in 
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Fig. 5. — The best fit values of the main pulse (PI), 
interpulse (P2) positions and distance between them 
in the phase plane versus the revolution number, in the 
20-100 keV energy band. The horizontal dashed lines 
show the one sigma (statistical only) confidence inter- 
vals for the averaged values. 

the radio timing ephemeris. 

4. Comparison with PCA 

For independent check of the hard X-ray results, 
we carried out the analogous analysis in the 2-20 keV 
energy band. We used the data of Crab observations 
with the PCA monitor quasisimultaneous with INTE- 
GRAL. We used the same radio ephemeris and the 
same pulse definition procedure. The phase positions 
of the X-ray main pulse relative to the radio one for 
79 PCA/RXTE observations are shown in the bottom 
panel of Figure [6] showing that also the main X-ray 
pulse leads the radio one. To determine the mean value 
of the time lag we approximated the PCA data with a 
constant and got 0.00933 ± 0.00016 in phase units or 
310 ± 6fis in time. Again we are quoting only statisti- 
cal errors but we keep in mind the systematic ~ 4-Q/is 
error that comes from the radio ephemeris uncertain- 
ties. Based on two measured time lags "X-ray/radio" 
and "hard X-ray/radio" we can conclude that the main 
X-ray pulse is leading the hard X-ray main pulse by 
35 + I6fis. This value differs only marginally from zero 



Fig. 6. — The main pulse maximum arrival phase in 
the 20- 100 keV energy band with SPI/INTEGRAL 
(1 point by time interval order of 100 ks), and in the 
2-20 keV energy band with PCA/RXTE (1 point — 
a few kiloseconds). The same radio ephemerides and 
fit procedure have been used. The horizontal dashed 
lines show the one sigma (statistical only) confidence 
intervals for the averaged values. 

(we provided the one sigma error) even taking into ac- 
count that in this case the radio error does not play a 
role since we used the same radio ephemerides for both 
measurements. Another quantity that is independent 
of any uncertainty in the radio timing ephemerides is 
the phase difference between the main pulse and in- 
terpulse. From our set of PCA/RXTE observations, 
we got a value of 0.40016 ± 0.00028 that is in a 
good agreement with the result obtained in hard X-rays 
(0.3989 + 0.0009, SPI/INTEGRAL). 

5. Radio delay evolution with energy 

We have measured accurately the delays of the main 
pulse arrival time in wide 2-20 keV (soft X-rays) and 
20 - 100 keV (hard X-rays) energy bands with respect 
to the radio main pulse arrival time. The 310 + 6 fis 
soft X-ray/radio delay derived in this paper based on 
the PCA data (Rots et al. 2004 provides even higher 
value 344 + 40 /us) is marginally higher than the radio 
delay with respect to 275 + 15 fis hard X-ray/radio de- 
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lay measured with SPI. Both values are also slightly 
higher of the radio/optical delay of 255 ± 21 fis de- 
rived from S-Cam optical observations (Oosterbroek 
et al. 2008) and the 241 ± 29 /us radio/gamma one 
(> 30 MeV, EGRET, Kuiper et al. 2003). To check 
whether this differences are real or not, we have made 
an additional analysis and investigated the behaviour 
of the radio delay with energy. The radio delay evolu- 
tion with energy have been observed in optical wave- 
length (Oosterbroek et al. 2008), that gave us an extra 
motivation. We have built folded curves in narrower 
energy bands. We split the 2-20 keV PCA energy 
band on three parts, while we used five energy chan- 
nels to cover the 20-300 keV energy band for SPI data. 
We also added the High-Energy X-Ray Timing Exper- 
iment (HEXTE/RXTE) data in four energy bands cov- 
ering the 20 - 250 keV energy range, allowing a direct 
crosschecking with SPI results. The position of the 
main peak has been determined as previously, from a 
fit with a composite model: a Gaussian function plus a 
constant background, in the phase interval 0.98 - 1 .01 . 
Figure |7]presents the evolution of the radio delay ver- 
sus energy for the RXTE and SPI data, together with 
the optical and y-rays points from Oosterbroek et al. 
(2008) and Kuiper et al. (2003). The decreasing trend 
of the radio delay with energy in the (2 - 300 keV) 
energy domain supports the reality of the delay mea- 
sured between the soft X- and hard X-ray main peaks, 
even though the individual error bars are large. When 
modelling this decrease by a simple linear law (dashed 
line in Figure|7]) we find that the radio delay decreases 
with a rate of ~ 0.6 + 0.2 fis/keV. The obtained Chi2 
of 7.1 for 10 dof compared to the value of 42.4 for 1 1 
dof for a constant model corresponds to a probability 
of 3.5 xlO -5 that its improvement is by chance. 

6. Discussion and conclusions 

We have investigated the pulse profile of the Crab 
pulsar between 2 and 300 keV with PCA/RXTE in- 
strument and SPI/INTEGRAL telescope. We found 
the strong indication that in this energy range the ra- 
dio delay is significantly decreasing with energy. The 
simplest explanation of such a behaviour is that the 
X - ray/y so f,-iay emission originates in a region ex- 
tended along the open magnetic field lines, with softer 
photons originating at higher altitudes while the time 
offsets represent simply the path-length differences. In 
this case, the time delay determines the characteristic 
azimuthal size of the emitting area as 54 + 18 km. 
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Fig. 7. — The radio delay in function of energy. The 
optical (S-Cam) and y-ray (EGRET) points are from 
Oosterbroek et al. 2008 and Kuiper et al. 2003, re- 
spectively. The data of the RXTE instruments and SPI 
are from this work. 

It is clear that the radio delay can not decrease lin- 
early through the full gamma rays range, since the 
EGRET point imposes a positive value. Moreover, 
contrary to we observe above 2 keV, the radio delay in- 
creases between optical and X-rays. It is interesting to 
note that the X-ray delay is consistent with the rate of 
-5.9+1 .9 us 1 1000 A derived from the 3920 - 8230 A 
optical waveband by Oosterbroek et al. (2008). 

Even thought the dependence of the main pulse 
phase position on energy is complex, we can ex- 
plain it with a rather simple scheme, where the pulsed 
emission consists of the superposition of two inde- 
pendent components having different phase distri- 
butions and energy spectra. Indeed, such an an- 
alytical model has been introduced by Massaro et 
al. (2000) to interprete the BeppoSAX data in the 
0.1-300 keV energy band (see also im provements 
Massaro et al. 2006t Campana et al. 20091) . In the op- 
tical up to y-ray domain, this model includes: an 
"optical", C and an "X-ray", C x . The fractional part 
of the C x component increases in the main pulse with 
increasing energy from 1 keV up to a 1 MeV, then de- 
creases up to ~ 10 MeV energy but is negligible below 
1 keV and above 10 MeV. C A /C ratio behaviour fol- 
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lowing the same law as the Bridge/Pi ratio presented 
in Kuiper et al. (2001). Thereby, in the X-ray band, 
the C x component shifts the maximum of the pulse I 
(C x + C ) emission toward the radio maximum on the 
phase plane. It explains that the radio delay decreases 
with energy in X-rays while it keeps identical values in 
the optical and y-ray wavebands where the C x com- 
ponent is negligible. Considering that the values of the 
radio delay are nearly the same in the optical wave- 
lengths and y-rays, where C a largely dominates over 
C x , we can suggest that the C component corresponds 
to a single emission mechanism and emission location 
from optical to y-rays. In this case, the ~ 250 fis radio 
delay indicates that the radio emission is produced in 
a region located closer from the Neutron Star by ~ 75 
km than the C production site. On the other hand, the 
C x component could be unrelated to the C one, with a 
different origin and/or source location. In this case, the 
source behaviour in X-rays would result from a super- 
position of (at least) two independent components and 
the measured values of the radio delay in this energy 
band would not have any direct physical explanation. 
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notes and discussion. We wish to thank to the INTE- 
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Table 2: The Crab pulsar ephemerides for the INTEGRAL observations listed in Table Q] We used the Crab Pulsar 



position given in Jodrell Bank: RA = 83.633217° and DEC = 22.014464° in the J2000 epoch 



Rev. 

N° 


T vaM (MJD), 
t'' u (MJD), 
$ p (sec) 


/ , Hz 
dfldt , 10 -10 sec" 2 
d 2 f/dt 2 , IQ- 21 sec" 3 


Atcfisec 
AY 


Rev. 

N" 


T valid (MJD), 
t' Q "< (MJD), 
t™ p (sec) 


/ - Hz 
df/dt , 10~ 10 sec" 2 
d 2 f/dt 2 , 10- 21 sec- 3 


Ato,fisec 
AT 


0043 


52671 -52699 
52685 
0.076659 


29.8092705147 
-3.7366060 
9.0 


18 
0.0005 


0300 


53431 -53461 
53444 
0.033023 


29.7847841837 
-3.7315793 
6.0 


-19 
-0.0006 


0044 


52671 -52699 
52685 
0.076659 


29.8092705147 
-3.7366060 
9.0 


-1 
-0.0000 


0352 


53584-53615 
53597 
0.029626 


29.7798524524 
-3.7299236 
9.0 


36 
0.0011 


0045 


52671 -52699 
52685 
0.076659 


29.8092705147 
-3.7366060 
9.0 


-69 
-0.0021 


0365 


53644 - 53675 
53658 
0.022656 


29.7778867428 
-3.7294045 
4.0 


-24 
-0.0007 


0102 


52852 - 52883 
52866 
0.018825 


29.8034282349 
-3.7350193 
-1.0 


8 

0.0002 


0422 


53796 - 53826 
53809 
0.005391 


29.7730221322 
-3.7278153 
4.0 


-19 
-0.0006 


0103 


52852 - 52883 
52866 
0.018825 


29.8034282349 
-3.7350193 
-1.0 


31.0 
-0.0009 


0464 


53917 -53948 
53931 
0.020193 


29.7690932051 
-3.7267461 
15.0 


4 

0.0001 


0124 


52913 -52944 
52927 
0.019997 


29.8014598690 
-3.7344076 
3.0 


-11 

-0.0003 


0483 


53971 -54009 
53993 
0.011812 


29.7670971393 
-3.7264807 
18.0 


162 
0.0048 


0125 


52913 -52944 
52927 
0.019997 


29.8014598690 
-3.7344076 
3.0 


-52 
-0.0016 


0541 


54160-54191 
54174 
0.020652 


29.7612711958 
-3.7245631 
5.0 


-25 
-0.0007 


0126 


52913 - 52944 
52927 
0.019997 


29.8014598690 
-3.7344076 
3.0 


-31 
-0.0009 


0605 


54344 - 54374 
54358 
0.010645 


29.7553516229 
-3.7225368 
13.0 


31 
0.0009 


0170 


53068 - 53074 
53071 
0.010486 


29.7969173441 
-3.7535080 
350.0 


-50.0 
-0.0015 


0665 


54526 - 54557 
54540 
0.008472 


29.7494993041 
-3.7206976 
13.0 


-46.0 
-0.0014 


0182 


53095 -53105 
53100 
0.008622 


29.7958809510 
-3.7414057 
200 


73 
0.0022 


0666 


54526 - 54557 
54540 
0.008472 


29.7494993041 
-3.7206976 
13.0 


-23 
-0.0007 


0184 


53105 -53115 
53110 
0.024904 


29.7955577865 
-3.7392771 
-100 


-22 
-0.0007 


0727 


54710-54741 
54724 
0.013181 


29.7435856030 
-3.7189687 
17.0 


54 
0.0016 


0239 


53255 - 53279 
53263 
0.021234 


29.7906210263 
-3.7333983 
7.0 


40 
0.0012 


0728 


54710-54741 
54724 
0.013181 


29.7435856030 
-3.7189687 
17.0 


69 
0.0021 


0247 


53279 -53311 
53293 
0.101225 


29.7896533932 
-3.7329692 
9.9 


-48.0 
-0.0014 


0774 


54863 - 54892 
54877 
0.000804 


29.7386704404 
-371735.72 
14.9 


-23 
-0.0007 
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